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Children born very preterm (VPT) are at risk for visual impairments, the main risk factors being retinopathy of 
prematurity and cerebral white matter injury, however these only partially account for visual impairments in 
VPT children. This study aimed to compare optic radiation microstructure and volume between VPT and term- 
born children, and to investigate associations between 1 ) perinatal variables and optic radiations; 2) optic radi- 
ations and visual function in VPT children. We hypothesized that optic radiation microstructure would be altered 
in VPT children, predicted by neonatal cerebral white matter abnormality and retinopathy of prematurity, and 
associated with visual impairments. 

142 VPT children and 32 controls underwent diffusion-weighted magnetic resonance imaging at 7 years of age. 
Optic radiations were delineated using constrained spherical deconvolution tractography. Tract volume and 
average diffusion tensor values for the whole optic radiations and three sub-regions were compared between 
the VPT and control groups, and correlated with perinatal variables and 7-year visual outcome data. 
Total tract volumes and average diffusion values were similar between VPT and control groups. On regional anal- 
ysis of the optic radiation, mean and radial diffusivity were higher within the middle sub-regions in VPT com- 
pared with control children. Neonatal white matter abnormalities and retinopathy of prematurity were 
associated with optic radiation diffusion values. Lower fractional anisotropy in the anterior sub-regions was 
associated with poor visual acuity and increased likelihood of other visual defects. 

This study presents evidence for microstructural alterations in the optic radiations of VPT children, which are 
largely predicted by white matter abnormality or severe retinopathy of prematurity, and may partially explain 
the higher rate of visual impairments in VPT children. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 



Abbreviations: AD, Axial diffusivity; BWSDS, Birth weight standard deviation score; CI, 
Confidence interval; CSD, Constrained spherical deconvolution; FA, Fractional anisotropy; 
GA, Gestational age; MRI, Magnetic resonance imaging; MD, Mean diffusivity; RD, Radial 
diffusivity; ROP, Retinopathy of prematurity; VPT, Very preterm. 
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1. Introduction 

Children born very preterm [VPT, <32 weeks' gestational age (GA)] 
are at high risk of visual impairments compared with children born at 
term (>37 weeks' GA) (Arpino et al., 2010), which range from relative- 
ly subtle complications such as reduced visual acuity to complete blind- 
ness (Mirabella et al., 2006). Complex visual perceptual processes that 
rely on the integration of visual information are also often affected 
(Geldof et al., 2012). Such visual impairments may contribute to learn- 
ing difficulties at school (Cooke et al., 2004). Major causes of visual im- 
pairment in VPT children include retinopathy of prematurity (ROP; 
pathologic vascular hyperproliferation in the immature retina (Palmer 
et al., 2005)) and overt cerebral white matter injury affecting the visual 
pathways (Ricci et al., 2006). However, visual impairments have been 
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reported even in VPT children who have no history of ROP or overt 
cerebral injury (Kozeis et al., 2012), indicating that other factors may 
be associated with visual impairment in VPT children. By studying the 
white matter within the optic radiations, elucidation of potential neuro- 
anatomical correlates for visual impairment in VPT children may be 
possible. 

One method that can be utilized to study the major tracts of the white 
matter is constrained spherical deconvolution (CSD) based tractography. 
CSD is an advanced diffusion-weighted magnetic resonance imaging 
(MRI) analysis method that is capable of modeling multiple fiber orien- 
tations, enabling robust delineation of white matter tracts (Tournier 
et al., 2007). Combined with the diffusion tensor model, non-invasive es- 
timates of the microstmcture and volume of specific white matter tracts 
can be made. Diffusion imaging has been applied to investigate the optic 
radiations in preterm infants in association with visual function assessed 
both in the neonatal period (Bassi et al., 2008; Berman et al., 2009; 
Groppo et al., 2012), and 6-20 months later (Glass et al., 2010). Howev- 
er, it is unknown whether further maturation and development of the 
optic radiations later into childhood strengthen or weaken relationships 
with visual function. The current study extends previous studies by in- 
vestigating the association between optic radiation microstmcture and 
volume and visual function at 7 years of age, when visual system struc- 
ture and function is approaching full development (Birch and 
O'Connor, 2001 ). Furthermore, it associates optic radiations with neona- 
tal MRI-detected white matter abnormality and ROP. It also extends pre- 
vious studies by including a term-born control group, having access to a 
large sample size, and utilizing advanced CSD based tractography. The 
specific aims of the current study were to: 1) compare the microstmc- 
ture and volume of the optic radiations between VPT children and con- 
trols at 7 years of age; 2) investigate relationships between perinatal 
variables and optic radiation structure in VPT children, particularly ROP 
and MRI-detected white matter abnormality; and 3) determine the rela- 
tionship between optic radiation structure in VPT children and visual 
function assessed contemporaneously. It was hypothesized that: 1 ) 
optic radiation structure would be altered in VPT children compared 
with controls; 2) neonatal white matter abnormality and ROP would 
be associated with altered optic radiation structure in VPT children; 
and 3) altered optic radiation structure would be associated with visual 
impairments in VPT children. 

2. Materials and methods 

2.3. Participants 

A prospective observational cohort of 227 children born <30 weeks' 
GA or < 1 250 g (which we will call the "VPT" group, acknowledging that 
it is different from the usual use of "VPT" as <32 weeks' GA) and 46 
healthy term controls (born > 37 weeks' GA) was recruited from the 
Royal Women's Hospital in Melbourne between July 2001 and Decem- 
ber 2003. 2 VPT infants died and 1 was subsequently found to have a 
congenital abnormality, leaving 224 VPT children eligible. Children 
with congenital abnormalities or chromosomal disorders were also ex- 
cluded. A total of 1 98 VPT children and 43 controls were followed-up at 
age 7 years. Of those who were followed-up, 160 VPT children and 36 
controls underwent MRI, but 22 of these either did not have full diffu- 
sion datasets acquired or scans were unusable due to movement 
artifact. Thus the final number of children in the current study was 
174 (64% of those recruited), including 142 VPT and 32 controls. 
Perinatal characteristics were similar between VPT participants 
who contributed optic radiation data and non-participants, with 
the exception of postnatal corticosteroid exposure (participants — 5.7%; 
non-participants —13.8%, p = 0.03). 

The study was approved by the Royal Children's Hospital Human 
Research Ethics Committee and was performed in accordance with the 
ethical standards laid down in the 1964 Declaration of Helsinki and its 
later amendments. Informed consent was obtained from all parents/ 



caregivers prior to inclusion in the study. Perinatal data were obtained 
by chart review at the time of hospital discharge. 

2.2. Magnetic resonance imaging 

At term-equivalent age, and T 2 weighted images were acquired on 
a 1.5 T scanner. Infants were immobilized in a vacuum fixation beanbag 
and scanned while sleeping, without sedation. A cerebral white matter 
abnormality score (range 0-17) was graded qualitatively, as previously 
described and validated (Kidokoro et al., 2013). 

At 7 years of age, MRI was performed without anesthesia or sedation 
on a 3 T scanner. weighted three-dimensional rapid gradient-echo 
images were acquired: TR, 1900 ms; TE, 2.27 ms; matrix, 256 x 256; 
FOV, 210 x 210 mm; 0.8 mm 3 isotropic voxels. Intracranial volumes 
were obtained from the T x weighted images using 'FreeSurfer' software 
(Buckner et al., 2004). Given that overall head size has been shown to 
influence white matter diffusion tensor values due to differential partial 
volume effects (Takao et al., 201 1 ), intracranial volumes were used to 
control for partial volume effects or other white matter microstructural 
differences related to head size. 

Two sets of echo-planar diffusion-weighted images were acquired; 
one with 25 non-collinear gradient directions and b-values ranging up 
to 1200 s/mm 2 (TR, 12000 ms; TE, 96 ms; matrix, 144 x 144; FOV, 
250 x 250 mm; 1.7 mm 3 isotropic voxels), and another with 45 gradient 
directions and a b-value of 3000 s/mm 2 (TR, 7400 ms; TE, 106 ms; ma- 
trix, 104 x 104; FOV, 240 x 240 mm; 2.3 mm 3 isotropic voxels). The 
b = 1200 s/mm 2 data were processed using 'ExploreDTF software 
(Leemans et al., 2009). Data were corrected for motion and eddy current 
induced distortions, incorporating re-orientation of the B-matrix. The dif- 
fusion tensor model was fitted, generating axial diffusivity (AD), radial 
diffusivity (RD), mean diffusivity (MD), and fractional anisotropy (FA) 
maps. CSD was applied to the b = 3000 s/mm 2 diffusion-weighted data 
using 'MRtrix' software (Tournier et al., 2012), creating a map of fiber ori- 
entation distributions in each voxel. A maximum harmonic order of 6 was 
used. 

2.3. Fiber tractography 

Probabilistic tractography was performed based on CSD's fiber ori- 
entation distributions. A seed region of interest was placed on an axial 
slice depicting the transition from the posterior limb of the internal cap- 
sule to the cerebral peduncle. The seed was positioned immediately 
lateral to the lateral geniculate nucleus within the white matter of the 
optic radiation at the apex of the arc around the lateral ventricles 
(Fig. 1 A), as described previously (Ciccarelli et al., 2003). A target region 
of interest was positioned on a coronal slice encompassing the entire 
white matter cross-section of the optic radiation, just anterior to the pri- 
mary visual cortex (Fig. IB) (Berman etal., 2009). A maximum fiber ori- 
entation distribution amplitude of 0.2 was specified, and tracts were 
constrained to pass in one direction from the seed to target region of in- 
terest (Fig. 1C). Binary tract images were produced by counting the 
number of streamlines per voxels, and discarding voxels containing 
< 6/1 000 streamlines to reduce the likelihood of including aberrant 
fibers in the tract volume. 

Diffusion metrics vary along the optic radiations due to variations in 
fiber density and geometry, partial volume effects and locations of inju- 
ry (Berman et al., 2009; Yeatman et al., 2012). Therefore, the tracts were 
divided into three equal sub-regions by cutting the tracts one-third and 
two-thirds of the way from their most anterior to most posterior points 
(Fig. ID, E and F). 

Diffusion values (FA, RD, AD, MD) within the whole optic radiations 
and sub-regions were obtained by averaging values from the co- 
registered b = 1200 s/mm 2 diffusion tensor maps within the binary 
tract images. Diffusion values were obtained from the co-registered 
b = 1200 s/mm 2 maps rather than directly from the b = 3000 s/mm 2 
maps because: 1 ) diffusion tensor values, particularly MD, vary according 
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Fig. 1. Regions of interest for optic radiation tractography. (A) Seed regions of interest were defined on an axial slice and positioned immediately lateral to the lateral geniculate nuclei 
within the white matter of the optic radiations at the apex of the arc around the lateral ventricles. (B) Target regions of interest were placed on a coronal slice, designed to encompass 
the entire white matter cross-section of the optic radiations as they passed into the primary visual cortices. (C) Optic radiation white matter tracts of a very preterm 7-year-old. 
(D), (E), (F) Locations of the anterior, middle and posterior sub-regions of the optic radiations respectively. 



to the strength and number of t>-values in the acquisition (Melhem et al., 
2000); 2) increasing b-values are accompanied by decreasing signal-to- 
noise ratios (Jones et al., 2012). As such, the b = 1200 s/mm 2 sequences 
may be expected to provide more accurate estimates of the diffusion ten- 
sor values than the b = 3000 s/mm 2 sequences. All co-registered 
images were visually examined for registration errors. 

Intra-class correlation coefficients for tractography repeated on 20 
randomly chosen participants were: tract volume = 0.88, FA = 0.93, 
AD = 0.82, RD = 0.98, MD = 0.95; indicating high intra-rater reliabil- 
ity by the single operator (D. Thai). 

2.4. Visual function examination 

Vision was assessed by a pediatrician masked to the clinical history 
of the children on the day of the 7-year MRI. Assessments were per- 
formed with the child's own (if any) corrective eyewear. Visual acuity 
was measured simultaneously for both eyes, using the standard Snellen 
eye chart, from a distance of 6 m. Results were recorded as Snellen dec- 
imals (Kniestedt and Stamper, 2003). Visual defects in either eye were 
recorded from previous opthalmological assessments and included stra- 
bismus, myopia, hypermetropia, cataracts, nystagmus, astigmatism, 
amblyopia or requirement for corrective lenses. Visual perception was 
assessed using the visual closure subtest from the Test of Visual Percep- 
tual Skills — Third Edition (Martin, 2006), which requires children to 
identify whole figures from incomplete images. Relationships between 
visual closure score and measures from both the right and left optic 
radiations of VPT children were assessed separately, as the right hemi- 
sphere is considered dominant for visual closure ( Wasserstein et al., 
2004). 

2.5. Statistical analysis 

Data were analyzed using 'SPSS v20' (Chicago, 1L). Exploratory anal- 
yses were initially performed to assess the distributions of optic 



radiation diffusion and volume data and visual outcomes for normality 
and outliers. 

Baseline characteristics of the VPT and control groups that contribut- 
ed 7-year optic radiation data were compared using t-tests, chi-squared 
tests or Mann-Whitney U tests, as appropriate. 

Data from the whole optic radiations and sub-regions were com- 
pared between the VPT and control groups, with and without adjusting 
for age at MRI and intracranial volume, using linear regression. 

Associations between clinically important perinatal variables 
and the optic radiation data of VPT children were investigated 
using multivariable linear regression: GA at birth; birth weight stan- 
dard deviation score (BWSDS) (Cole et al., 1998); MRI-detected white 
matter abnormality score; female sex; severe ROP (maximum stage in 
either eye > 3); postnatal corticosteroid exposure. Each perinatal 
variable was firstly entered into separate regression models. Second- 
ly, all perinatal variables were entered simultaneously into the re- 
gression models, along with age at MRI and intracranial volume, to 
investigate the effect of each perinatal variable independent of the 
effects of the other perinatal variables, age at MRI and intracranial 
volume. 

Within the VPT group, relationships were investigated between: 
1 ) The optic radiation variables and visual acuity score using linear re- 
gression; 2) The optic radiation variables and the likelihood of visual de- 
fect of any kind in either eye using logistic regression; 3) The optic 
radiation variables and visual perception score using linear regression. 
These analyses were performed with and without adjusting for age at 
MRI, intracranial volume and the previously defined perinatal variables. 
Given that visual perception may be related to overall cognitive func- 
tioning, the third stage of analysis also excluded children with an intel- 
lectual impairment [full-scale IO_ score < 70, as determined using the 
Wechsler Abbreviated Scale of Intelligence (WASI) (Wechsler, 1999)]. 

Given the multiple comparisons, the results were interpreted by 
considering overall patterns and magnitudes of differences, rather 
than by focusing on individual p-values. 
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Table 1 

Perinatal characteristics and 7-year visual outcomes of the very preterm ( VPT) and control cohorts. 
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Mean difference (95%C1) 


p 


Age at imaging (years), mean (SD) 


7.5 (0.2) 


7.6 (0.2) 


-0.08 (-0.17, 0.02) 


0.11 


Full-scale 1Q. score, mean (SD) 


98.8 (12.8) 


109.6(11.5) 


-10.8 (-15.6, -5.9) 


< 0.001 


Total intracranial volume (cm 3 ), mean (SD) 


13313 (121.8 ) c 


1421.8(104.8) 


-90.5 (-136.5, -44.5) 


< 0.001 


Total white matter volume (cm 3 ), mean (SD) 


403.6 (46.8) e 


444.6 (42.9) f 


-41.1 (-59.6, -22.6) 


< 0.001 


Visual acuity, mean (SD) 


0.99 (0.21 ) s 


1.12 (0.22) f 
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0.002 


Visual closure total raw score, mean (SD) 


5.2(3.1) c 
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Odds ratio (95% CI) 


p 


Visual defect of any kind in either eye, n (%) 


33 (23.6)' 


4 (12.5) 


2.2 (0.7, 6.6) 


0.17 


Strabismus, n (%) 


0 C01 


0 (01 






Astigmatism, n (%) 


9 (6.3) 


2 (6.3) 


1.0(0.2,4.9) 


0.99 


Myopia, n (%) 


8 (5.6) 


2 (6.3) 


0.9 (0.2, 4.4) 


0.89 


Amblyopia, n (%) 


1 (0.7) 


1(3.1) 


0.2 (0.01, 3.6) 


0.25 


Hypermetropia, n (%) 


5 (3.5) 


1 (3.1) 


1.1 (0.1, 10.0) 


0.91 


Cataract, n (%) 


0(0) 


0(0) 






Nystagmus, n (%) 


0(0) 


0(0) 






Corrective lenses prescribed, n (%) 


9 (6.3) 


1(3.1) 


2.1 (0.3, 17.2) 


0.48 


BWSDS = birth weight standard deviation score, CI = confidence interval, GA 


= gestational age, ROP = 


retinopathy of prematurity, SD 


= standard deviation. * Indicates not possible to 



calculate as there is a zero value in one or more cells. 
a BWSDS computed relative to the British Growth Reference data (Cole et aL, 1998). 

b Postnatal dexamethasone, usual dose 0.15 mg/kg per day for 3 days, reducing over 10 days: total dose 0.89 mg/kg. 

c n = 141. 

d n = 123. 

e n = 124. 

r n = 30. 

g n = 133. 

h n = 31. 

' n = 140. 



3. Results 

Table 1 reports sample characteristics of the 142 VPT children and 32 
controls with optic radiation data at age 7 years. Age at the 7-year MRI 
and the proportion of males was similar between the VPT and control 
groups, but many other perinatal characteristics differed as expected. Of 
note, intracranial volume was lower in the VPT children, and the VPT 
group exhibited poorer visual acuity and visual perception (visual closure). 

3.3. Very preterm children vs. controls 

There was higher diffusivity within the optic radiations of the VPT chil- 
dren than controls, particularly in the middle sub-region, even after 
adjusting for age at MRI and intracranial volume (Fig. 2A). There was 
evidence that optic radiation volume was smaller in the VPT children 
compared with the controls in the posterior sub-region in an unadjusted 
analysis [mean difference (95% confidence interval (CI)) = —0.4 (—0.8, 
— 0.04) cm 3 , p = 0.03], although this group difference was less marked 
after adjustment for age at MRI and intracranial volume (Fig. 2B). 

3.2. Perinatal predictors 

Within VPT children, higher neonatal white matter abnormality 
scores were associated with higher AD, RD and MD in the optic radia- 
tions, particularly in the middle sub-region. Higher neonatal white mat- 
ter abnormality scores were also associated with lower optic radiation FA 



in VPT children, both overall and in the anterior and middle sub-regions 
(Fig. 3 ). These associations remained after adjusting for age at MRI, intra- 
cranial volume and the remaining perinatal variables (data not shown). 
There was also evidence that severe ROP was associated with higher 
AD, RD and MD in VPT children, particularly in the middle sub-region 
(Fig. 3); adjustment for age at MRI, intracranial volume and the remain- 
ing perinatal variables had little effect (data not shown). There was little 
evidence of associations between GA at birth, BWSDS, sex or postnatal 
corticosteroid exposure and optic radiation diffusion values in VPT chil- 
dren (Fig. 3). 

Female sex was associated with lower optic radiation volume in VPT 
children, particularly in the posterior sub-region (Fig. 3E); adjustment 
for age at MRI, intracranial volume and the remaining perinatal vari- 
ables had little effect (data not shown). There was little evidence of as- 
sociations between GA at birth, BWSDS, postnatal corticosteroid 
exposure, neonatal white matter abnormality or severe ROP and optic 
radiation volume in VPT children (Fig. 3E). 

3.3. Visual outcome 

Within the VPT group there was little evidence of associations 
between optic radiation diffusion values and visual acuity, except higher 
FA in the anterior sub-region was associated with higher (better) visual 
acuity scores. This association remained after adjusting for age at MRI, 
intracranial volume and the previously defined perinatal variables 
(Fig. 4A). Increased optic radiation volume was associated with better 
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Fig. 2. Mean differences in optic radiation variables [(A) diffusion values within the optic radiations, (B) optic radiation volume] in very preterm (VPT) compared with term-born control 
children. Results presented are adjusted for age at MR1 and intracranial volume. AD = axial diffusivity (x 1CP 3 mrf/s), CI = confidence interval, FA = fractional anisotropy, MD = - 
mean diffusivity (xlCT 3 mm 2 /s), RD = radial diffusivity (xlO" 3 mm 2 /s). 



visual acuity for all regions, with similar results after adjusting for age at 
MRl, intracranial volume and perinatal variables (Fig. 4B). 

Reduced optic radiation FA both overall and within the anterior and 
middle sub-regions, of VPT children was associated with increased odds 
of visual defect in either eye, a finding that remained after adjustment 
for age at MRI, intracranial volume and perinatal variables (Fig. 5A). 
There was also some suggestion that increased RD was associated 
with increased odds of a visual defect. There was little evidence of asso- 
ciations between optic radiation volume and visual defect in VPT chil- 
dren (Fig. 5B). 

There was some suggestion that lower RD and MD, particularly in 
the anterior sub-region, of both the left and right optic radiations of 
VPT children was associated with higher (better) visual closure scores, 
with some attenuation of these associations after adjusting for age at 
MRI, intracranial volume and perinatal variables (Fig. 6A and B). There 
was little evidence of associations between right or left optic radiation 
volume and visual closure in VPT children (Fig. 6C). These findings 
were similar when the only child with a full-scale IQ score <70 was 
excluded (data not shown). 

4. Discussion 

In summary, optic radiation diffusivity, particularly in the middle 
sub-region, was higher in the VPT children than controls. Neonatal 
white matter abnormality and severe ROP were associated with higher 
diffusivity within the optic radiations of VPT children, particularly with- 
in the middle sub-region. Female sex was also associated with lower 
optic radiation volume in VPT children. Furthermore, in VPT children 
increasing optic radiation FA, particularly in the anterior sub-region, 
was associated with better visual acuity and decreased likelihood of 
other visual defects. 

Diffusivity was higher in the optic radiations in the VPT group com- 
pared with controls. This finding suggests that the microstructure of 
the optic radiations is altered in VPT children compared with controls. 
Such changes in VPT children may reflect changes in white matter fiber 
diameter or density, myelination or membrane permeability (Jones 
et al., 2012). The findings may also suggest that optic radiation develop- 
ment is altered in VPT children compared with term-born children, 
which is supported by previous studies reporting that in typically devel- 
oping children AD, RD and MD decrease over time in the white matter 
(Mukherjee and McKinstry, 2006). The detection of diffusion differences 



between VPT children and controls in the middle sub-region to a greater 
extent than in the anterior and posterior sub-regions may relate to 
variations in axonal architecture along the tracts, or localization of 
injuries at specific points along the tracts (Berman et al., 2009; 
Yeatman et al., 2012). Additionally, there may be reduced influence 
by confounding factors such as crossing fibers in the middle sub- 
region, giving greater sensitivity to find group-wise differences. 
Optic radiation volumes, particularly in the posterior sub-region, 
appeared smaller in the VPT children. However, this finding was 
largely explained by reduced intracranial volume and/or confound- 
ing from age at MRI. This finding is consistent with Gimenez et al., 
who found little difference in optic radiation volume between VPT 
and control adolescents (Gimenez et al., 2006). 

Within the VPT group, higher MRI-detected neonatal white 
matter abnormality scores predicted increased diffusivity and re- 
duced FA in the whole optic radiations, and particularly the mid- 
dle sub-region. This suggests that optic radiation microstructure 
is particularly altered in VPT children who had higher white mat- 
ter abnormality scores in the newborn period. Similarly, there 
have been previous reports indicating that early white matter ab- 
normalities correlate strongly with altered white matter diffusion 
tensor values in preterm children and adolescents, including 
lower FA and higher AD, RD and MD (Bonifacio et al., 2010; 
Feldmanet al., 2012). 

VPT children with severe ROP in infancy displayed higher diffusivity, 
particularly in the middle sub-region, than VPT children with no or mild 
ROP (<stage 3), suggesting an association between severe ROP and 
altered optic radiation microstructure. It is possible that visual depriva- 
tion from the ROP-affected eye may directly affect visual pathway 
development. Supporting this, reduced optic nerve fiber thickness was 
recently found in preterm children with severe ROP compared with 
term controls (Akerblom et al., 2012). Relationships between severe 
ROP and altered optic radiation microstructure may also reflect associa- 
tion with a common etiological factor, such as hypoxic-ischemic insult, 
which has been implicated in the causal pathways to both ROP (Stout 
and Stout, 2003) and cerebral white matter injury (Volpe, 2009). In con- 
trast to our findings, Bassi et al. reported similar optic radiation diffusion 
tensor values between preterm infants with and without ROP at term 
equivalent age (Bassi et al., 2008). However, Bassi et al.'s study included 
only 8 infants with mild ROP (stage 1 or 2 ) that had resolved by the time 
of MRI. In the current study, 30 VPT children had severe ROP ( > stage 3 ) 
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Fig. 3. Relationships between perinatal variables and optic radiation variables in the VPT children [(A) fractional anisotropy, (B) axial diffusivity, (C) radial diffusivity, (D) mean diffusivity 
and (E) tract volume]. Results presented are regression coefficients, representing the difference in optic radiation per unit change in the predictor, from an unadjusted analysis. AD = axial 
diffusivity (x 10~ 3 mm 2 /s), BWSDS = birth weight standard deviation score, CI = confidence interval, FA = fractional anisotropy, GA = gestational age at birth (weeks), MD = mean 
diffusivity (x 10~ 3 mm 2 /s), PCS = postnatal corticosteroid exposure, RD = radial diffusivity (x 10~ 3 mm 2 /s), ROP = severe retinopathy of prematurity, WMA = white matter abnor- 
mality score, (J = regression coefficient. 



in infancy, enabling detection of microstructural differences between 
VPT children who had no or mild ROP and those with severe ROP in 
infancy. 

Female sex was associated with lower optic radiation volume in VPT 
children, particularly in the posterior sub-region, even after adjustment 
for intracranial volume, age at MR1 and the remaining tested perinatal 
variables. A possible explanation is that relative white matter volume 
increases more rapidly in males than females throughout childhood 
and adolescence, particularly in the frontal and occipital regions 
(Lenroot et al., 2007). 

In the VPT children, increasing FA in the anterior optic radiation sub- 
region was associated with better visual acuity and decreased likelihood 
of other visual defects. Similarly, previous work has reported correla- 
tions between diffusion values within the optic radiations of preterm 
neonates and functional visual outcome (Bassi et al., 2008; Berman 
et al., 2009; Glass et al., 2010). Together, these findings may suggest 
that structure-function relationships exist between optic radiation dif- 
fusion values and visual functions relying on visual cortical processing. 
The current study suggests that relationships between optic radiation 
structure and visual function in VPT children persist beyond infancy, 
and up to 7 years of age. 



There was evidence of associations between RD and MD in both the 
right and left optic radiations of VPT children and our measure of visual 
perception (visual closure). However these associations may have been 
explained by age at MRI, intracranial volume and/or the previously de- 
fined perinatal variables. Additionally, similar patterns of associations 
were observed for the right and left optic radiations, despite reports 
that the right visual cortex is dominant for visual closure ( Wasserstein 
et al., 2004). Right and left optic radiation diffusion values are highly 
correlated, therefore an association between microstructure and visual 
closure in one hemisphere is likely to be present in the other. 

A limitation of this study was that a proportion of children in the 
cohort did not have useable diffusion-weighted images. Furthermore, 
the participants in this analysis tended to be healthier, with a smaller 
proportion of VPT participants exposed to postnatal corticosteroids 
compared with non-participants, and our findings may have been 
stronger if the entire cohort had been able to be assessed with MRI at 
7 years of age. 

There are limitations associated with using Snellen charts to mea- 
sure visual acuity. Particularly in patients with low vision, the accuracy 
of Snellen acuity recordings may be affected by chart features such as 
non-uniform changes in letter sizes from line to line, variations in the 
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Fig. 4. Relationships between optic radiation variables [(A) diffusion values within the optic radiations, (B) optic radiation volume] and visual acuity score in VPT children. Solid lines 
represent regression coefficients from an unadjusted analysis and broken lines represent regression coefficients adjusted for age at MR1, intracranial volume and clinically important peri- 
natal variables (gestational age at birth, birth weight standard deviation score, sex, white matter abnormality score, severe retinopathy of prematurity and postnatal corticosteroid expo- 
sure). AD = axial diffusivity (xl0~ 3 mm 2 /s), CI = confidence interval, FA = fractional anisotropy, MD = mean diffusivity (xlfj~ 3 mm 2 /s), RD = radial diffusivity (xl0~ 3 mm 2 /s). 
P = regression coefficient. 



numbers of letters per line, and inconsistencies in the spacing between 
letters and lines. These and other disadvantages of Snellen charts have 
been discussed in detail previously (Kaiser, 2009; Kniestedt and 
Stamper, 2003). The accuracy with which visual acuity is measured 
may be improved in future studies by using, for example, Early Treat- 
ment Diabetic Retinopathy Study (ETDRS) charts rather than Snellen 
charts. Future studies may also benefit from assessing a wider range of 
visual functions, including stereopsis, contrast sensitivity and visual 
fields, in order to provide a more holistic view of visual system function- 
ing in VPT children and the potential role of the optic radiations. 



There are also limitations in using diffusion tensor measures to study 
the white matter. The tensor model enables detection of changes in sev- 
eral microstructural properties of the white matter (e.g. fiber density, 
myelination, etc.), but cannot sensitively distinguish between such fac- 
tors. Moreover, diffusion tensor values may be influenced by the intra- 
voxel orientational coherence of fibers and by partial volume effects 
from non-white matter tissues (Jones et al., 2012). Recent work has 
demonstrated that FA (Alexander et al., 2001) and MD (Vos et al., 
2012) are lower in crossing than single fiber regions, and that RD and 
AD may also be affected (Wheeler-Kingshott and Cercignani, 2009). 
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Fig. 5. Relationships between optic radiation variables [(A) diffusion values within the optic radiations, (B) optic radiation volume] and risk of visual defect in VPT children. Solid lines rep- 
resent regression coefficients from an unadjusted analysis and broken lines represent regression coefficients adjusted for age at MRI, intracranial volume and clinically important perinatal 
variables (gestational age at birth, birth weight standard deviation score, sex, white matter abnormality score, severe retinopathy of prematurity and postnatal corticosteroid exposure). 
AD = axial diffusivity (xl0~ 3 mm 2 /s), CI = confidence interval, FA = fractional anisotropy, MD = mean diffusivity (xl0~ 3 mm 2 /s), RD = radial diffusivity (xlCT 3 mm 2 /s), 
P = regression coefficient. 
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Fig. 6. Relationships between optic radiation variables [(A) diffusion values within the left optic radiation, (B) diffusion values within the right optic radiation, (C) optic radiation volume], and 
visual closure scores in VPT children. Solid lines represent regression coefficients from an unadjusted analysis and broken lines represent regression coefficients adjusted for age at MRI, in- 
tracranial volume and clinically important perinatal variables (gestational age at birth, birth weight standard deviation score, sex, white matter abnormality score, severe retinopathy of pre- 
maturity and postnatal corticosteroid exposure). AD = axial diffusivity (x 10~ 3 mm 2 /s), CI = confidence interval, FA = fractional anisotropy, MD = mean diffusivity (x 1CT 3 mm 2 /s), 
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Future work may use the recently developed 'Apparent Fibre Density' 
(Raffelt et al., 2012), which is based on CSD's fiber orientation distribu- 
tions, to improve interpretation of diffusion changes in crossing fiber 
regions. 

Tract thickness, orientation and curvature may also affect the diffu- 
sion tensor measures, in part by influencing the extent of partial volume 
effects (Vos et al., 201 1 ). The current study attempted to control for par- 
tial volume effects by adjusting for intracranial volume in statistical 
analyses. Along with intracranial volume, other measures and methods 
have been proposed to control for partial volume effects, such as total 
white matter volume or regional tract volume, and there is some uncer- 
tainty as to which is most appropriate (Metzler-Baddeley et al., 2012). 

Additionally, despite using CSD ductography to obtain robust tract 
delineations (Tournier et al., 2007), this study was unable to consistent- 
ly track Meyer's loop. 

5. Conclusions 

VPT birth is associated with altered optic radiation microstructure in 
childhood, particularly within the middle sub-region, relative to birth at 
term. Such microstructural alterations in VPT children may be predicted 
by complications resulting from VPT birth such as MRl-detected white 
matter abnormalities and severe ROP. Moreover, altered optic radiation 
microstructure in VPT children, particularly in the anterior sub-region, 
is associated with reduced visual acuity and increased likelihood of 
other visual defects. These findings may partially explain the high 
rates of visual impairments in VPT children. By further elucidating the 
neural correlates of visual impairments in VPT children, the current 
findings may ultimately aid in improving the functional visual outcome 
of VPT children. 
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